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ABSTRACT 

We present a two-dimensional kinematic analysis for a sample of simulated binary 
disc merger remnants with mass ratios 1:1 and 3:1. For the progenitor discs we used 
pure stellar models as well as models with 10% of their mass in gas. A multitude 
of phenomena also observed in real galaxies are found in the simulations. These 
include misaligned rotation, embedded discs, gas rings, counter-rotating cores and 
kinematic misaligned discs. Using the 2D maps we illustrate projection effects and 
the change in properties of a merger remnant when gas is included in the merger. 
We find that kinematic peculiar subsystems are preferably formed in equal mass 
mergers. Equal-mass collisionless remnants can show almost no rotation, regular 
rotation or strong kinematic misalignment. The inclusion of gas makes the remnants 
appear more round(l:l) and axisymmctric(3:l). Counter-Rotating Cores (CRCs) are 
almost exclusively formed in equal-mass mergers with a dissipational component. 
3:1 remnants show a much more regular structure. We quantify these properties 
by applying the kinemetric methods recently developed by Krajnovic et al. This 
work will help to understand observations of elliptical galaxies with integral field 
spectrographs, like SAURON. 

Key words: methods: analytical - methods: N-body simulations - galaxies: elliptical 
and lenticular, cD - galaxies: formation - galaxies: evolution - galaxies: fundamental 
parameters - galaxies: kinematics and dynamics 



1 INTRODUCTION 

N-body simulations of mergers of disk galaxies have been 
very succesful in reproducing many observational fea- 
tur es of low and intermediate luminosity elliptical galax - 
ies (|Toomre fc Toomre|[l973 : lBarnes|[l993lHernquistll 19921 ) . 
while the most massive ellipticals are probably a re- 
sult of collisionless early-type mergers or mixed mergers 
ijKhoehfar fc Burkert|[20o3 : lNaab et al.ll2006l ). However, the 
observational information obtained from galaxies has been 
greatly enhanced with the advent of Integral Field Uni t 
(IFU) instruments such as SAURON (|Bacon et al.ll2001^ ■ 
OASIS and CMOS. While mostly long-slit data was used to 
constrain the kinematic properties of early-type galaxies un- 
til recently, 2-dimensional data contains more information, 
and makes e.g. the detection of a second distinct kinematic 
component, which might be aligned with an angle not cov- 
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ered by a slit in a conventional observation, much easier or 
even possible at all. 

The SAURON collaboration presented the first rep- 
resentative sample of 2 -dimensional observat ions of El- 
liptical and SO galaxies l|Emsellem et al.l |2004| . henceforth 
EM04). They presented 2D maps of line-of-sight velocities, 
dispersions and the third and fourth Gauss-Hermite mo- 
ments. These maps showed a multitude of features such 
as kinematically decoupled components (KDCs, see also 
iDavies et al.ll200ll ). central velocity dispersion dips, counter- 
rotating disks and many others. Dynamical modelling of 
the kinematic data revealed that galaxies with significant 
rotation arc more anisotropic than slow-rotating ellipticals 
(Cap pcUari et al. 2005). An interesting new correlation was 
found that the slow-rotating ellipticals host large KDCs 
which have stellar populations of the same age as the main 
body, in contrast to fast-rotating ellipticals which have much 
smaller KDCs and are younger on average (jMcDermid et all 
[2OO61I . This correlation must contain some important infor- 
mation on the formation history of early-type galaxies. 

Numerical studies of merging galaxies in the past 
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showed that many pecuhar kinematical features occur nat- 
uraUy in merger remnants, especially when a dissipational 
component was present in the progenitor galaxies. The dis- 
sipative component which can cool falls to the center and 
changes the morpholgy of the remnant drastically by re - 
arranging its orbital content ([Barnes fc HernauistI Il996l ). 
Disk like counter-rotating cores can also be formed from 
infalling gas, which has a spin vecto r with the opposite 
sign than the main rotating body (|Hernauist Sz Barnes! 
I1991I ). however also dissipationless scenarios have been pro- 
posed for counter-rotating stellar populat ions in early-type 
galaxies ( Balcclls & Gonzalez 1998:Balcc lls fc QuinnI [l990l : 
iBendo fc Barnes! (2000!). IBarnes! (|2002l ) showed that the 
gaseous material which does not fall into the center can 
form large scale gas discs which can be warped or form bars. 
Polar-ring like features have been succesfuUy produced by 
several authors e ither in binary mergers jB ckki 199^), tidal 
accretion event s llBournaud Sz Cornbesll2003! ) or in a cosmo- 
logical context l|Macci6 et al.l 2006! ). 

The aim of this work will be to compare 2-dimensional 
observations of a representative sample of 1:1 and 3:1 disc- 
disc merger remnants to the observations of the SAURON 
galaxy sample (EM0 4). The only simi la r stud y of this type 
was carried out by iBendo fc Barnes! (|2000! ). who identi- 
fied peculiar 2-dimensional kinematic features in collision- 
less merger models, such as as an orthogonally decoupled 
core, counter-rotating populations and misaligned rotation. 
Their sample contained equal-mass mergers as well as merg- 
ers with mass ratios of 3:1 and they found that equal-mass 
mergers produce a larger variety of kinematical features. But 
as pointed out before, coUisionless mergers are unlikely to 
be the only formation channel for early-type galaxies and 
it is unclear how the gas will influence 2D maps of such 
remnants. Gas infall also will vary from merging geometry 
to mergin g geometry a s the graviational torques on the gas 
will vary i|Barnes!l200'3 '). 

W e will examine the merger sample of iNaab et al.l 
l|2006! ). henceforth NJB06, which includes for every merg- 
ing geometry a coUisionless merger and a merger where each 
progenitor galaxy had 10% of its disc mass in gas. In total 
we study a sample of 96 mergers which should hopefully give 
us a similar variety of kinematic features than found in the 
SAURON sample. As we can compare the mock obsevations 
of each coUisionless remnant with its gaseous counterpart, 
we can easily assess which kinematical feature was already 
present in the dissipationless remnant and which was caused 
by the influence of the gas.. 

Although 2-dimensional features can many times be in- 
tuitively grasped by visual inspection of the maps, this is 
hardly a quantifiable way to compare to simulations. Ke- 
centlv lKrainovic et al.! l|2006! ). henceforth K06, developed a 
method to describe the properties of 2-dimensional maps 
of arbitrary moments of the line-of-sight-velocity distribu- 
tion(LOSVD). They termed this method kinemetry which 
works much in analogy to the common photometric fitting 
of isophotes. With this method one can determine, for ex- 
ample, in an easy way position angle, amplitude and shape 
of the rotation of a galaxy, or merger remnant, and iden- 
tify kinematic subsystems more easily. The shape of the 2- 
dimensional maps are of course viewing angle dependant, 
and we try to take this into account by analysing the pro- 
jections along the three principle axes. However, we are not 



carrying out a comprehensive statistical study in this work. 
We will use kinemetry throughout this work to quantify 
the kinematic features which arise through merging of disc 
galaxies. 

In Section[2]we introduce the simulation sample. We ex- 
plain the observational method to obtain 2D maps of various 
moments of the LOSVD in Section [3] The maps and their 
kinemetric analysis are presented in Section [4] We compare 
these results to 2D-kinematics of observed galaxies in Sec- 
tion [5] and summarize our findings in Section [5] Velocity 
maps and kinemetry of the whole sample are shown in the 
Appendix. 



2 SIMULATIONS 

The coUisionless simul ations are a subse t of th e simulations 
discussed in detail by iNaab fc Burkerj (|2003l ) , henceforth 
NB03. The simulations with gas, which we use for our anal- 
ysis are identical to the ones presented in NJB06. In the 
following we give only the most important simulation pa- 
rameters: 

The progenitor disc galaxies were constructed in 
dynamical equilibrium using the method described by 
Hcrnquist (|l993l ). The system of units was: gravitational 
constant G=l, exponential scale length of the larger disc 
in the merger hd = 1 (the scale height was hz — 0.2) and 
mass of the larger disc Ma ~ 1. The discs were exponential 
with an additional spherical, non-rotating bulge with mass 
Mt — 1/3, a Hernquist density profile (Hcrnquist 1990) and 
a scale length Vb = 0.2h, and a pseudo-isothermal halo with 
a mass Md = 5.8, cut-off radius rc — lOh and coUision- 
less core radius 7 — Ih. The parameters for the individual 
components were the same as for the coUisionless mergers 
presented in NB03. For this study we have re-simulated the 
full set of 1:1 and 3:1 mergers with an additional gas com- 
ponent in the disc. We replaced 10% of the stellar disc by 
gas with the same scale length and an initial scale height of 
hz.aas = 0-lhz. This is i n agreement with recent results from 
iKhochfar fc^fl^ | 2006! ) in which they show that merging 
progenitor galaxies at low redshifts contain on average 10% 
gas for ellipticals more massive than 10^" solar masses. The 
gas was represented by SPH particles assuming an isother- 
mal equation of state, P = c^p, with a fixed sound speed of 
Cs = 0.039 in velocity units corresponding to Ca ~ 10km/ s if 
scaled to a Milky Way type galaxy. The N-body/SPH sim- 
ulations were performed using the hybrid N-body/SPH tree 
code VINE (Wetzstein et al. in prep.) with individual time 
steps. 

The galaxies approached each other on nearly parabolic 
orbits, in agreement with prediction s from cosmological sim- 
ulations (IKhochfar fc Burkert!l2006l ). with an initial separa- 
tion of 30 length units and a pericenter distance of 2 length 
units. The inclinations of the two discs relative to the or- 
bit plane were ii and 12 with arguments of pericenter wi 
and LU2- In selecting unbiased initial parameters for the disc 
inclinations we followed the procedure described bv lBarne^ 
(1998). The initial orientations for the discs were the same 
as in NB03, Table 1. The merger remnants were allowed 
to settle into dynamical equilibrium for approximately 10 
dynamical timescales after the merger was complete. Then 
their equilibrium state was analyzed. 
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Table 1. Properties of the modeled merger remnants. Special 
kinematical features observed and initial orientations of the pro- 
genitor discs are shown. 
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The total sample comprises of 96 merger remnants of 
which half of them were run with a gaseous component. We 
adopt the following nomenclature for the mergers presented: 
The first two numbers give the mass ratio of the progenitor 
discs, i.e. 11 or 31. The letter 'C stands for analysing a col- 
lisionless merger (disc and bulge particles), 'S' stabds for 
analysing only the stellar component of a gaseous merger 
(disc and bulge particles) and 'GS' analysing the total lu- 
minous part of a gaseous merger (disc, bulge and gas parti- 
cles). The last number denotes the merging symmetry. For 
example, the abbreviation 11C5 means: a collisionless equal 
mass merger with merging symmetry 5. The mergers we have 
analysed in detail, a comment regarding their kinematic pe- 
culiarity and the detailed orbit parameters of the merger 
they resulted from, can be found in Table [1] 



3 2-DIMENSIONAL DATA ANALYSIS 
3.1 LOSVDs 

Every remnant was projected along the long axis (X- 
axis, YZ- projection), the intermediate axis (Y-axis, XZ- 
projection) and short axis (Z-ajds, XY-projection) of the 
moment of inertia tensor defined by the 40% most tightly 
bound stellar particles. For the 2D analysis we binned parti- 
cles within the central 3 length units on a grid of 48 x 48 cells. 
This corresponds t ypically to 2-3 effectiv e radii depending 
on projection (see Naab &: Truiilld |2006| for the exact de- 
termination of Te). To include seeing effects we created for 
every luminous particle 10 x 10 pseudo-particles with identi- 
cal velocities on a regular grid with a total size of 0.125 unit 
lengths centered on the original particle position. The mass 
of the original particle was then distributed to the pseudo- 
particles weighted by a Gaussian with a standard deviation 
of 0.1625 unit lengths. Thereafter the pseudo-particles were 
binned on a 48 x 48 grid. 

For the kinematic analysis we binned (mass weighted) 
all pseudo-particles falling within each grid cell in veloc- 
ity along the line-of-sight. The width of the velocity bins 
was set to a value of 0.1 for line-of-sight velocities vios in 
the range —4 ^ wios ^ 4. This resulted in 80 velocity bins 



over the whole velocity interval. Using the binned velocity 
data we constructed line-of-sight velocity profiles (LOSVD) 
for each bin of the 2D grid. Subsequently we parameter- 
ized deviations from the Gaussian shape of the velocity 
profile using Gauss-Hermite basis functions (|Gerhardlll993l : 
Ivan der Marel fc Franx|[l993l ). The kinematic parameters of 
each profile (aat, vm, H3, H4) were then determiiied si- 
multaneo usly by least squares fitting (jCretton et al.ll200ll : 
iNaab fc Burkertll2001bD . 



3.2 Kinemetry 

Surface brightness and line-of-sight velocity are both mo- 
ments of the stellar distribution function. The surface 
brightness is its zeroth-order and the line-of-sight velocity 
the first-order moment. However, there is a fundamental 
difference between the 2D fields of both moments: as the 
surface brightness is an even moment its iso-contours are 
point-symmetric (and therefore closed) and as the line-of- 
sight-velocity is an odd moment its iso-contours are point- 
anti-symmetric and therefore open. While isophotal shapes 
were found to be very close to perfect ellipses, such a choice 
is not obvious for iso- velocity contours. K06 found empir- 
ically that the kinematics of early-type galaxies resemble 
very closely those of inclined discs, which in many cases can 
be well described by circular motion. If you would project a 
circular orbit which moves in the plane of the galaxy on the 
sky, then it would follow a simple cosine law: 



V{R,Tp) ^Vo + Vc{R) sin i cost/;, 



(1) 



where R is the semi-major axis length of the projected circle 
on the sky, Vb the systemic velocity, Vc is the ring circular 
velocity, i is the ring inclination and ip is the azimuthal angle 
measured from the major axis of the ellipse (their Eq. 5). 

However, as the authors state, elliptical galaxies are 
spheroids and not inclined discs and therefore their assump- 
tions will be violated in nature. The next step is to see how 
far galaxies deviate from the simple cosine law. This can 
be achieved by the means of a Fourier series. A map of a 
moment of the LOSVD can be generalized to 

JV 

K{a,ilj) = Ao(a)-l-^^„(a)sin(n?/>) + B„(a)cos(ni/'), (2) 

n = l 

where t/j is the eccentric anomaly and a the semi-major axis 
of the ellipse (their Eq. 6) . K06 find that only few terms are 
needed (up to third order) to find a good description of the 
2D data. 

In the case of an even moment the fitting procedure is 
similar to what is done in conventional photometry, while in 
case of an odd moment their algorithm tries to minimize all 
terms of the harmonic expansion except Bi (a) cos('!/)), which 
according to their Ansatz carries the major contribution to 
the map of an odd moment of the LOSVD. In the first stage 
of the fit, these terms are minimized on a grid of position 
angles (F) and flattenings [q = 1 — e). After the ellipse, 
which describes the velocity profile best, has been found, 
the rotation curve is extracted again and Fourier- analysed 
to arbritrary order. 

We follow the authors in using the following four pa- 
rameters to quantitatively describe kinematic 2D maps: The 
already mentioned F, which measures the alignment of the 
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rotation and q which can be interpreted as the opening angle 
of the iso-velocity contours. They argue that although the 
amplitudes of the separate cosine and sine terms in principle 
describe different properties of the maps, it is better for sys- 
tems which are not axisymmetric to collect the amplitudes 
of sine and cosine terms of the same order (their Eq. 10): 

fc„ = y^Al + Bl. (3) 

For n=l we get the term fci, which represents the amount 
of line-of-sight bulk rotation in case of a velocity map. The 
authors found that the ratio k^/ki, where fcs is the first 
higher-order term of the expansion which has not been fit- 
ted to the map and expresses the deviations from regular 
rotation, is sensitive to the presence of multiple kinematic 
components in a galaxy. 

For a detailed description of their analysis method we 
refer the reader to the original paper K06. 



4 RESULTS 

We show in this section the most important 2-dimensional 
maps of our sample in order to cover a broad range of differ- 
ing kinematic features. The selection of these maps is again 
guided by, but not restricted to, the results of EM04. We will 
use the results of our previous study on the orbital structure 
of merger remnants l|jesseit et al.ll2005l . henceforth JNB05) 
and a kinemetric analysis for the interpretation of the maps. 
The apparent long axis of the remnants will be horizontally 
oriented for all the maps presented in this paper. In gen- 
eral we show maps for the collisionless stellar particles (C), 
the stars of the simulations with gas (S) and the stars in 
combination with gas (GS), as an indication of what the 
star would look like in the case of late star formation. The 
newly 'formed' stars would not have the same stellar mass- 
to-light ratio than the 'old' stellar population. To account 
for this effect we adopt an ad hoc mass-to-light ratio follow- 
ing common stellar population models and modify the maps 
accordingly. During the merger a significant fraction of the 
gas falls to the center (see NJB06) and overpowers the signal 
in the central bins and limits in some cases our interpreta- 
tion of the results at very small radii. All maps are shown 
in the original SAURON colour scheme as devised by Cap- 
pellari & Emsellem (2001) and included in the Kinemetry 
distribution. 



4.1 Surface Density Maps 

It has been shown by a number of investigations that the 
presence of gas dur ing a merger event makes the remnant 
more axisymmetric l|Barnes fc Hernquistlll996l ) . N JB06 have 
shown that this effect is stronger for 1:1 remnants which are 
intrinsically more triaxial (or prolate) than 3:1 remnants. 
Fig.[T]shows this effect for the 2D stellar surface density dis- 
tribution of a collisionless 1:1 remnant and its counterpart, 
simulated with gas. Just by eye it is obvious that this par- 
ticular remnant becomes much rounder if gas is included. Its 
intrinsic b/a changes from 0.83 to 0.95, while c/a increases 
even more significantly from 0.51 to 0.81. This is one of the 
few remnants in which the major axis tube fraction increases 
compared to the collisionless case from 5% to 11%. Next to 



minor axis tubes, also outer major axis tubes support po- 
tentials of spherical shape. So an increase of this fraction 
is to be expected (for a detailed description of orbit classi- 
fication in merger remnants see JNB05). The effect of gas 
on the surface density maps of 3:1 remnants is in general 
hardly detectable (see Fig. ^ as the collisionless remnant is 
already rather axisymmetric. 

4.2 Rotation of Collisionless Remnants 

In Fig. Owe show the line-of-sight- velocity field of two colli- 
sionless 1:1 remnants and one collisionless 3:1 remnant 31C6 
with regular rotation. The properties of 3:1 remnants do not 
change significantly with the initial disc orientations of the 
progenitors. Therefore this example is typical for this mass 
ratio. The remnant 11C3 shows also regular rotation, despite 
of being an equal-mass merger, albeit with a lower ampli- 
tude than 31C6. 11C3 is the 1:1 remnant with the highest 
minor axis tube fraction, but it is rather an exception, as in 
general the rotational structure of equal-mass merger rem- 
nants is more complicated and depends on the initial disc 
geometry. In the most extreme case a 1:1 merger remnant 
can exhibit almost no rotation within in an effective radius, 
like 11C5, which has the highest box orbit fraction of all 
merger remnants. 

The kinemetry shows (Fig.|3l bottom row), as indicated 
before, that 11C3 and 31C6 have a rising bulk rotation, i.e. 
fci curve, while 11C5 has almost no amplitude in the center. 
We extract the box orbit population in each remnant, to see 
how this changes the maps. Mainly we observe now particles 
which have been classified as minor axis tubes, which have 
a significant amount of L^. In every case the amount of 
rotation increases (bottom row, right panel) and also the iso- 
velocity contours appear to be more closed. Still, all three 
maps are quite different in appearance. That means that 
in the merging process also the shape of the z-tubes is re- 
arranged. They are considerately more 'puffed up' in equal- 
mass mergers than in un-equal mass mergers. 

We want to note that the 2D velocity maps of 1:1 merg- 
ers are sometimes too complex to be fitted adequately with 
ellipses. Very low rotation causes the ellipse fitting algorithm 
to extract erratic position angles. We therefore restrict our- 
selves ot constant position angle of 0° and constant elliptic- 
ity of g = 0.5. See also discussion in Sec. 14. 4|) 

4.3 The Influence of Gas on vios and /13 

In Fig. |4] we show the 2D line-of-sight-velocity map of 
a 3:1 remnant and the corresponding distribution of /13. 
For the collisionless remnant (left panels) vi^b and /13 are 
correlated inside 0.5 re and anti-correlated in the outer 
parts. This behaviour is typical for collisionle ss remnants 
(IBendo fc Barnes|[2OO0l : lNaab fc Burkertl2001al : NJB06) and 
is not consistent w ith observations of elliptical galaxies 
(|Bender et al.lll994^ . The correlation of v and /is at the cen- 
ter is caused by the superposition of box orbits and minor 
axis tubes in collisionless merger remnants. This happens 
because the peak of the total velocity profile at a radius 
close to the center is significantly moved to lower velocities, 
with respect to the short axis tube population, because box 
orbits, having zero mean angular momentum, peak normally 
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Figure 1. 2D surface density maps for the equal mass remnant with geometry 8. The projections along the three principal axes for the 
collisionless remnant {11C8, top row) are compared to the stellar component of the corresponding remnant with gas (11S8, bottom row). 
11C8 is clearly triaxial whereas 11S8 is nearly round due to the influence of the addit ional gas component on the distribution of stars. 
One unit of lengths corresponds to roughly one effective radius llNaab fc Truiilldl2006l) . Colours indicate surface density in logarithmic 
scale. 




-1 D 1 -1 a 

Figure 2. Same as Fig. [T]for the 3:1 merger with geometry 6. There is no obvious influence of gas on the stellar distribution. 



at zero mean velocity. Therefore the normally steep leading 
wing becomes broader (see NJB06 for a more detailed dis- 
cussion). When we remove the box orbit component and 
observe this modified remnant, the correlation at the center 
has disappeared (Fig. U top middle panel). 

If gas is present during the merger the correlation at 
the center disappears even if only the stellar component 
is considered (|4j top row, right panel), leading to an over- 
all anti-correlation which is even stronger if we include the 



gas component in the analysis (two bottom rows). The 
gas influences the kinematic structure of the remnant two- 
fold: first, it is suppressin g the population of box orbits 
(jSarnes fc Hernauistlll99^ . because the gradient of the po- 
tential in the center is very steep and orbits which come 
close to the center, i.e. b ox-like, are scattered on c entropho- 
bic, i.e. tube-like orbits (jCerhard fc Binnevl[l985l ). Second, 
the gas particles themselves settle into a disk-like configura- 
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Figure 3. 2-dimensional line-of-sight velocity field for the coUisionless 1:1 remnant 11C5 (top row), the 1:1 remnant 11C3 (second 
row) and the 3:1 remnant 31C6 (third row). The bulk rotation is indicated in the bottom row. 31C6 (green triangles) and 11C3 (black 
circles) show regular rotation, while 11C5 (red upside down triangles) has a large box orbit component and rotates very little. How the 
2-dimensional maps of the remnants would look like with the box orbits removed is shown in the last column. The rotation curve without 
box orbits is plotted with thick symbols and for comparison the original rotation curve with thinner symbols (bottom right). The net 
rotation increases significantly when the box orbits are removed. 



tion. Both processes work towards a, v — correlation more 
in agreement with what is found in real galaxies (NJB06). 

As the hght carried by particles on short axis tubes is 
important for the correct correlation, and as the gas parti- 
cles move mostly in the XY-plane we want to test how a 
different stellar mass-to-light-ratio would affect our results. 
We endorse a very simple method by giving more weight to 
the gas particles relative to the stellar particles. To assign 
realist ic M/L ratios we use the results o f iBruzual fc CharlotI 
for a stellar population of solar metallicity Z=0.02, 
Salpeter IMF and in the K-Band. A 2 Gy old population 
of such stars would have a M/Lk, which is 2.527 lower and 
5 Gy old population a M/Lk which is 1.336 times lower 
than a 10 Gy old population. The TO Gy population' being 
the case when stars and gas particles have the M/Lk, i-e. 
a normal 'GS' observation. When we multiply the masses 



with the correction factors we see in Fig. |3] (two bottom 
rows) that the effect on the maps is minimal, but visible. 
The peak velocity in the 2 Gy observation is higher and the 
iso- velocity contours somewhat more closed, but our conclu- 
sions are unchanged. 

4.4 Kinematic Twists and Misalignment 

A galaxy is kinematically misaligned if its rotation axis is 
not aligned with its photometric major axis. Because of the 
strong kinematic misalignment the assumption of a constant 
position angle is of course wrong. To solve this problem we 
follow the alternative procedure, proposed by K06, and ex- 
tract the rotation curve by averaging over circles (see Ap- 
pendix A of K06). Despite this simplification the authors 
point out that fci and fcs/fci are still useful indicators for 
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Figure 4. 2d maps of velocity and hi for the 3:1 remnant with geometry 19. For the coUisionless remnant (31C19) uios and /13 is 
correlated in the inner parts and anti-correlated in the outer parts (top left). This correlation vanishes, when the box orbits are removed 
(top middle). Gas similarly changes the orbital setup and no correlation is visible in the stellar component of the gaseous remnant (top 
right). If the gas is 'ad hoc' turned into stars the anti-correlation becomes even morevisible (two bottom rows). Accounting for the lower 
M/L of the gas particles increases the peak velocity but changes hardly the structure of the maps. 



the bulk motion and deviations from regular motion. But 
because a circle has no defined position angle, we have to 
recover the phase angle 



bi = arctan 



(4) 



where the phase angle 0i represents the direction of the bulk 
rotation (their Eq. 10). 

By convention signifies a (fi ~ 180° misaligned rotation 
and — 90° alignment with the major photometric axis. 
As shown in the top row of Fig. [S]the coUisionless remnant 
11C12 is at the center almost maximally misaligned with 
a <j>i = 180° and slowly changes to about (pi = 150° at 



larger radii, i.e. the remnant has a kinematic twist (KT). 
Tliis reflects the fact that the fraction of minor axis tubes 
increases with increasing radius, which skews the rotation 
more towards the major photometric axis. Nearly the same 
happens for the gaseous remnant (top row middle panel), 
only that the strongly misaligned region is smaller, as minor 
axis tubes are already populated at smaller radii. The KT 
is slightly more pronounced, to about <f>i = 125°, between 
an outer disc-like component and a population of major axis 
tubes in the center. 

In the second row of Fig. [5] we show an example of an- 
other coUisionless remnant, IICIO, with slightly lower mis- 
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Figure 5. Velocity maps of two 1:1 merger remnants with geometry 10 and 12 of the stars of the collisionless and gaseous remnants 
(column one and two). Kinematic position angle is shown in the third column for collisionless (black circles) and gaseous (red triangles) 
remnant. Misaligned rotation (180°) is indicated by a horizontal dashed red line and aligned rotation by a continuos black line. The 
maps of the collisionless remnants show strong kinematic misalignment (left column), which diminishes at larger radii. The inclusion of 
gas affects both remnants differently. While the center of IISIO exhibits a strong kinematic twist and then almost aligned rotation is the 
kinematic position angle of 11S12 hardly changed (last column). 



aligned rotation ranging from 01 — 160° to 125°. Inter- 
estingly, if re-simulated with gas (IISIO), the stellar rem- 
nant shows a changed kinematic structure with respect to 
the collisionless remnant. The misalignment is almost com- 
pletely removed. It is evident that the gas has to overcome 
a stronger centrifugal barrier in the remnant 11S12 which 
partially prevents it to sink to the center and reorganize the 
orbital structure of the remnant. In remnant IISIO this is 
not true and the gas manages to accumulate in the XY plane 
re-arranging the shape of the remnant. This is an indication 
that very strong misalignment can survive some degree of 
dissipation while small kinematical misalignment is wiped 
out more easily. 

4.5 The Flattening of Iso- Velocity Contours 

K06 argued with the example of four early-type galaxies 
taken from the SAURON sample that different kinematic 
flattenings, q, show to what extent a galaxy is an axisym- 
metric rotator, in which case q remains constant with radius 
(see their Fig. 7) . Two of their galaxies show peculiar flat- 
tenings, one with a radially rising q and one with very wide 
opening angles in the center (NGC 2549), but dropping to 
q = 0.4 at larger radii. 

As the q profile of NGC 2549 reminds us of flattenings 
found in 3:1 remnants, we will study in this section the influ- 
ence of gas on the q-parameter in our 3:1 merger remnants. 
Their kinematics are much simpler than in 1:1 remnants, 
because it is much more difficult to populate major axis 
tube orbits in a 3:1 merger remnant and hence their kine- 
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Figure 6. Flattening, q, of the iso-velocity contours vs radius. 
The q parameter has been averaged at each radial bin for all 
3:1 remnants for four cases: 3:1 collisionless merger XZ and YZ 
projection and 3:1 gaseous merger XZ and YZ-projcction. Error 
bars indicate the rms at each bin. For comparison we plot the 
kinemetric fit of a progenitor disc galaxy 



matic twists are almost negligible. Also, all of them show 
significant rotation in most of the projections, which is a 
precondition to apply the kinemetric analysis. We analyze 
the YZ and the XZ projection, i.e. edge on projections, of 
every remnant. To get an overview of how the g-parameter 
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changes with projection and with the inclusion of gas, we 
calculate the mean q for every radial bin and compare four 
subsamples: 3:1 coUisionless YZ and XZ-projection and 3:1 
dissipational mergers YZ and XZ-projection. We can see in 
Fig. |6]that the iso-velocity contours of the XZ-projection of 
the coUisionless remnants have in general a smaller open- 
ing angle in the outer regions than for the YZ-projection, 
reflecting the more extended, boxy shape of the minor axis 
tubes for the YZ projection (NJB06). The stellar remnants 
whose progenitors had a gaseous component, are almost ax- 
isymmetric and therefore the flattenings are not very dif- 
ferent between the two projections. It is, however, evident 
from the plot that the q of the stars in the gaseous rem- 
nants is significantly lowered at radii smaller than two effec- 
tive radii. Interestingly for large radii we find slightly larger 
opening angles in the gaseous remnants as compared to the 
XZ-projection of the coUisionless remnant. This is probably 
because at large radii the surviving part of the progenitor 
disc dominates the rotation and the flattening of the initial 
disc is very low as indicated in Fig. [6] (blue line). 

4.6 Systems with Polar Rings 

A particularly interesting 1:1 remnant is shown in Fig. [T] 
In this case the gas forms a large scale gas ring which is 
inclined by F = 95° with respect to the photometric major 
axis of the galaxy and a central gaseous disc. The kinemetric 
analysis reveals the different amount of rotation in the gas 
and in the stars. As the gas ring is very narrow we extract the 
rotation curve along an ellipse with a, q — 0.3 and with the 
same alignment, i.e. F = 95°, as the ring. As apparent from 
the bottom left panel of Fig. [7] already in the coUisionless 
remnant does a significant fraction of stars move on polar 
orbits. The gas also seems to have a higher line-of-sight bulk 
motion than the stellar component. 

The polar ring is the result of the initial conditions of 
this particular merger for which initially the discs of the 
progenitors are inclined to each other by nearly 90°. The 
properties of this remnant are similar to those of polar ring 
galaxies. Bokki (1998) has shown that such systems can form 
in galaxy merging events. Our simulations differ insofar that 
both progenitor discs have gas and that we ran a larger 
sample of disc inclinations. Only in two out of 16 1:1 mergers 
was a significant amount of gas deposited in a polar ring. 
Typically this is just 10% of the initial gas mass or 1% of the 
total baryonic mass of the system. The gas in such a polar 
ring stems from one of the merging partners only. We do not 
find any significant polar ring structures in 3:1 mergers. 

Observed polar ring galaxies, though, have a compa- 
rable amount of mass in the polar ring and in the central 
component, typically an SO galaxy (see e.g. ISparke fc Coxl 
I2OO0I ). Our simulations already start with too little gas (10% 
of the total mass) to form such a syst em. Alternatively, s uch 
systems can form by cold accretion (jMaccio et al.|[2006f ) in 
a cosmological context. 

4.7 Counter- Rotating Cores 

The detection of Kinematically Decoupled Components 
(KDCs) at the center of early-type galaxies is an important 
indication for the hierarchical assembly of elliptical galax- 
ies. A subset of these are central stellar sub-systems, which 



rotate in the opposite direction to the outer stellar body, 
so called Counter- Rotating Cores (CRCs). However, the ex- 
act process of how material, supposedly a left over from a 
merger, retains memory of its original an g ular m omentum 
vector is a matter of debate. iKormen d3 1 19841 ) proposed 
that a dwarf satellite under certain pre-conditions can sur- 
vive the tidal field of the more massive progenitor and sink 
to the center. iBalcells fc Quin3 ([19901) then showed in N- 
body simulations that such a formation scenario for a CRC 
is viable f or certain encounter geome tries. It was demon- 
strated by iHernquist fc BarnesI (Il99ll ) that dissipation can 
play an important role i n forming central counter - rotating 
components (but see also lBalcells fc Gonzalejl 19981 ). iBarned 
(2002) elaborated on these findings, but pointed out that 
for the same encounter geometry (ii = 71°,ciJi = 90°, 22 ~ 
— 109°,a;2 = 90°), but for a wider orbit no CRC is formed. 
The exact reason remains elusive. However, kinematically 
misaligned gaseous components were commonplace in their 
simulations, especially for retrograde encounters. 

Similarly we find kinematically misaligned gas, counter- 
rotating or otherwise, in about 50% of the 1:1 merger rem- 
nants, but in only 1 out of 32 3:1 remnants. Our best exam- 
ple for a CRC in the gas is the merger remnant 11GS16, 
which originates from a very similar merging geometry 
(ii = 60°,a)i = 90°, i2 = -109 °, 0^2 = 71°) than the 
original IHernquist fc BarnesI (|l99ll ) calculation. In Fig. [5] 
we show edge on 2-dimensional velocity maps of this rem- 
nant. Kinemetry is uniquely suited to analyse maps with 
KDCs. K06 showed a very instructive example of a set of 
two-component test models in which they vary the align- 
ment of the central component from aligned (AF = 0°) to 
counter-rotating (AF = 180°) (their Fig. 2) and proofed 
that the size and alignment of the KDC can be easily mea- 
sured by kinemetry. We take advantage of this and extract 
the kinemetry of the stellar component of the gaseous rem- 
nant, of the gas amd stars with identical M/L and of gas 
and stars with the gas particles having a lower M/L ('2Gy 
old population'). 

Observing the Gas-CRC at later or earlier times is cru- 
cial for the detectability in the maps. While kinemetry does 
pick up the signal in the TOGy' observation (Fig|8] middle 
column), the '2Gy' map shows that the possibility for an 
observer to find this CRC would be much better. Interest- 
ingly the old stellar component is marginally dragged along 
into counter-rotation, although a weak signal is detectable. 
This CRC is relatively small, about 10% of the effective ra- 
dius and the outer disk is lop-sided. The CRC is formed in 
this particular remnant due to a preservation of the orig- 
inal spin alingment of the progenitor discs. The rotation 
of particles which belonged to either disc 1 or disc 2 in the 
corresponding, i.e. originating from the same merging geom- 
etry, coUisionless remnant are counter-rotating. Apparently 
in the dissipational remnant does the gas component of one 
progenitor galaxy lose more angular momnentum than the 
other and settle in a fast spinning disk in the center. 

However, the most clearly visible and largest CRCs are 
found in the old stellar component of equal-mass merger 
remnants. 11S6 is such an example and its maps and kineme- 
try is shown in Fig. [9] We can neither detect a CRC in the 
coUisionless remnant (top left), nor in the gas (top right), 
but only in the 'old stars' (top middle). The CRC has a 
size of approximately half an effective radius and shows all 
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Figure 7. Top row: Edge on projection (along the intermediate axis) of the 2D Velocity fields of 11C9, 11S9 and 11GS9. The rotation 
curve was extracted along ellipses with q = 0.3 and F = 95° and are overplotted in the top middle panel.. Gas settles in a ring with 
a diameter of about two effective radii, inclined by 95° with respect to the main stellar body. The signature is similar to polar ring 
galaxies. Very little mass is contained in the polar gas ring itself. Bottom row: The polar rotational structure is already present in the 
collisionless remnant (bottom left panel, the bulk rotation of the corresponding map above is marked by thicker symbols). The gas ring 
shows a higher rotation than the stellar component in the gaseous remnant. 



tell-tale signs in the kinemetry, e.g. 180° change of kine- 
matic position angle. It is also visible in the showing 
an anti-correlation with vias as would be expected from a 
axisymmetric system with regular rotation. To study this 
phenomenon closer we examined the four distinct stellar 
components from which the particles which build up the 
CRC can consist: disc of galaxy 1, bulge of galaxy 1, disk 
of galaxy 2 and bulge of galaxy 2. In Fig llOl (top) we illus- 
trate that the most strongly counter-rotating components 
are the particles from the disc and the bulge belonging orig- 
nally to the same galaxy (the particles of galaxy 2 lost most 
of their rotation and lie between those two extremes). We 
found it instructive to examine the evolution of the angle 
between the spin vectors of the bulge particles of galaxy 1 
and the other three stellar components, repectively the two 
other gas components. We calculate the angular momentum 
vector of each component with respect to its own center 
of mass. In the same Fig. 1111 left panel, we see that the 
alignment of the spin of the counter-rotating bulge changes 
suddenly around t = 90, which is shortly after the second, 
and final, encouter (around t = 70). The bulge seems to 
be impulsively removed from its own disk. Later evolution 
shows that counter-rotation is slightly increasing and the 



angle difference stabilizes at —150° to —180°. Only a sig- 
nificant part of the second bulge is at later times dragged 
into counter-rotation as well. This process is not occuring 
in all the remannts. The spin vector of all components, gas 
or stellar, of the remnant 11GS3 (Fig. 1111 right panel) are 
aligned and no CRC is visible in the maps. 

The formation of a CRC in such a context opens up 
questions if this is a robust result if realistic modes of star 
formation are included or if the subsequent merging history 
preserves such a KDC. Some aspects of forming a CRC in 
the old stellar component are favourable in light of recent 
observations and are discussed in Sec. \S\ These issues are 
beyond the scope of this paper and we try to address them 
in a subsequent paper. 

4.8 Line-of-sight velocity dispersion 

The 2-dimensional velocity dispersion maps of elliptical 
galaxies of the SAURON sample show a richness of features, 
which are even less understood than velocity maps in terms 
of their formation history. Some of the galaxies show low dis- 
persion regions, mostly in the center. Drops in a are indica- 
tive that dissipation has played a role in the formation of the 
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Figure 8. Top row. 2-dimensional velocity maps of the 11GS16 witii a counter-rotating gas dislc in the center. The old stellar component 
left, as well as the gas and stars with the same M/L middle and lower M/L for the gas particles right are shown. The kinemetry is 
extracted along ellipses with q = 0.3 and F = 0°. Second row: The peak of the rotation curve (fci) and the visibility of the CRC is 
improving from left to right. The CRC is hardly visible in the old stellar component. Third row: The 180° kinematic position angle 
change is best visible in the 2Gy old population map. A second hump is visible at r 0.4 and is a sign for the lopsidedness of the outer 
gas disk. Bottom row: Again the 2Gy map shows the clearest transition signature in kc,/k\. 



galaxy. It is indeed a very good indicator because lower ve- 
locity dispersion means that locally entropy has been taken 
out of the system and that would be hard to imagine in a 
pure coUisionles process. However, colhsionless mergers can 
still produce non-trivial velocity dispersion maps as is shown 
in Fig. 1121 The merger remnant llCl has a high cr disk-like 
feature, which is a result of the peculiar merging geometry 
(ii =Q°,LO\ = 180°, 12 — 0°,aj2 ~ 0°), i.e. both progenitors 
are in the same plane and have anti-aligned angular momen- 
tum. The two superposed counter-rotating 'disks' in the fi- 
nal remnant maximize the line-of-sight velocity dispersion. 
In contrast the merger 11C3 shows a more typical flat dis- 



persion profile. Both remnants have a dispersion dip in the 
center, but this signature is a lef t-over from the c old centers 
of the initial Hernquist bulges (|Hernauistl Il990) . Also the 
velocity maps are quite different for both remnants. While 
the two counter-rotating components in llCl compensate 
each other, resulting in almost zero vios, does the remnant 
11C3 show normal rotation. The high g disk phenomenon 
is probably rare, as llCl is the only remnant in our sample 
who has this feature, which is no surprise as the merging 
set-up needs to be rather fine-tuned. 

Merger remnants with gas show a much wider variety 
in their velocity dispersion maps. We illustrate this with 
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Figure 9. Top row: Projected 2D line-of-sight velocity in tlie XZ projection for tlie collisionless remnant (11C6), the stellar remnant 
of tlie counterpart with gas (11S6) and including the gas (11GS6) component assuming that all gas turned into stars after the merger 
was complete (from left to right). In the middle panel a counter-rotating stellar core (CRC) is clearly visible. In the right panel the 
kinematic signature of the thin disc appears. The ellipses, which have a. q = 0.5 and F = 0°, along which we extracted the rotation 
curves are overplotted. Second Row: Rotation curves, fci, extracted along an ellipse with q = 0.5 and F = 0° centered on the photometric 
major axis. The curves are shown for the collisionless remnant (11C6), gas and stars of the dissipative remnant (11GS6) and the stellar 
component of the gaseous remnant (11S6), which has a Counter-Rotating Core (CRC). Third Row: For 11S6, the ratio k^/ki shows a 
strong peak at the transition radius between the CRC and the outer spheroid, as found in the test examples of K06. The other curves 
show no distinct features. Fourth Row: Reconstructed angular direction 4>\ of the rotation, k\. 11S6 shows a CRC with almost exactly 
180° misalignment in its rotation (i.e. counter-rotation). Bottom row: Two dimensional distribution of ha for the same remnants. Strong 
signatures for the counter-rotating core and the thin disc (middle and right panel). 
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Figure 11. Left: Evolution of tlie angle between the angular momentum vector of bulgcl particles of remnant 11S6 and all other 
components: Dl(diskl), B2(bulge2), D2(dik2), Gl(gas diskl) and G2(gas disk2). Alignment is indicated by the thick horizontal black 
line at 0°. The spin of the Bulgel is suddenly changed shortly after the second and final passage and stays at almost —180° with respect 
to all other components. Only the particles of bulge 2 seem to evolve secularly and become more counter-rotating with time. Right: The 
same properties plotted for remnant 11GS3, which has no CRC. The spin angles of all dynamical components are aligned shortly after 
the merging process. 
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Figure 12. Top row: Velocity dispersion of two coUisionless 1:1 merger remnants llCl and 11C3. The remnant llCl has a high velocity 
dispersion disk-like component, as can be seen from the kinemetric profiles (black circles, right panel), significantly higher than 11C3 
(green triangles). Bottom row: Because llCl originated from an in-plane merger of two disk with opposite spins is almost no net rotation 
visible. 11C3 shows regular rotation. 



the example of a 1:1 and a 3:1 merger remnant (Fig. I13|) . 
The dissipation in the 3:1 remnant (bottom row) leads to 
a cold disk-like component, which imprints a dumbbell like 
velocity dispersion structure. In general gaseous disks in un- 
equal mass mergers settle into the equatorial plane, extend 
to larger radii and have more rotation than in equal-mass 
mergers. The kinemetry shows that line-of-sight dispersion 
is lower than in the coUisionless remnant at almost all radii 
(same figure, bottom right). Quite the contrary happens in 
the merger 11GS4, where a fraction of the gas particles ac- 
cumulate into a high dispersion ring at about an effective ra- 



dius from the center. Because the biggest fraction of the gas, 
like in all simulations, settles into a compact, dynamically 
cold, central disc, the dispersion map shows a double-peak 
structure (Fig. 1131 top row). The peaks being edge-on cuts 
through the torus-like structure of the gas at this radius. 
The dissipation at this radius actually causes an increase 
the velocity dispersion with respect to the coUisionless rem- 
nant (same figure, top right). 

Also face-on (XY) views of gaseous remnants have dis- 
tinctive features, such as ring-like depressions in the velocity 
dispersion (top row of Fig. I14|) . The most interesting aspect 
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Figure 13. Top row: The velocity dispersion maps of tlie XZ projection for the 1:1 merger with merging geometry 4. The velocity 
dispersion of the coUisionless merger, 11C4, shows a slight depression of aios at the center. The dispersion field of the stellar component 
of the gaseous remnant, 11S6, shows an off-axis double peak, which is even stronger if the gas particles are taken into account (11GS4). 
The kinemetry velocity dispersion extracted along an ellipse with q = 0.3 is shown in the most right column for 11C4 (black circles), 
11S4 (green triangles) and 11GS4 (red upside down triangles). Inclusion of the Bottom row: In the 31GS19 the gas settles partly into a 
cold stellar disc, which causes a V-shaped depression in the velocity dispersion. 
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Figure 10. Line-of-sight velocity averaged in a thin slit along 
the major axis at the end of the simulation of remnant 11S6. 
The particles originally belonging to the disk of galaxy 1 (black 
filled circles) show a clear counter-rotation with respect to the 
particles orignally belonging to the bulge of the same galaxy (open 
diamonds). 



is that this feature is also imprinted by the gas in the old stel- 
lar component. The low dispersion is not more pronounced If 
gas and star particles are observed simulataneously (11GS2), 
however the central gas component appears even colder in 
the map. The depression in a corresponds to a positive /14 



indicating a LOSVD that is more peaked than a Gaussian 
(same figure, bottom row). 



5 COMPARISON TO OBSERVATIONS 

The previous sections highlighted the diversity of 2- 
dimensional maps of various moments of the LOSVD. We 
want to discuss in this section, how these results relate to 
observations of real early-type galaxies, e.g. to the results of 
EM04. Unfortunately this comparison is mainly qualitative, 
although we will loosely refer to the kinemetric categories 
as suggested by K06. For each category we will refer to the 
galaxies which are the principal examples for this category 
and comment on their possible merger origin. 
Low-Rotation: The kinemetry of the remnant 11C5 (Fig. 

shows (scaled) rotation of below 10 km/s inside one ef- 
fective radius, which is extremely low-rotation. There are 
many galaxies in the SAURON sample that are not rotat- 
ing outsi de a fast rotating cent ral KDC, the so called slow 
rotators (|Cappellari et alj|2005l ). We find three examples of 
galaxies with a rotation below 20 km/s: NGC 4374, NGC 
4486 and NGC 5846. All these objects are round and have 
elliptical or slightly boxy isophotes. Their high central ve- 
locity dispersion indicates that they are fairly massive and 
that a merger of two very massive late-type galaxies would 
be needed to destroy the rotation in the progenitors and 
increase the mass accordingly. The scenario of merging el- 
fiptical galaxies would be even more favourable , as already 
the p rogenitors would have very little rotation (|Naab et al.l 
l2006h . 
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Figure 14. Projected 2D velocity dispersion for remnant 11S2 (top row) the tlie 2D distribution of /14 for the same projections (bottom 
row). A depression of a at about 0.5re is clearly visible. In the same region /i4 changes sign and becomes positive indicating a LOSVD 
more peaked than a Gaussian. 



Kinematic Twists and Misalignment: A KT is indi- 
cated by a smooth change of the kinematic position angle. 
We find KTs only in 1:1 mergers which are violent enough 
to scatter stars from minor axis tubes onto major axis tube 
orbits. More than half of the 1:1 merger remnants show, at 
least at the center, iso-velocity contours twisting more than 
30°. The frequency and amplitude in the simulated equal- 
mass merger remnants seems to be higher than for the KTs 
found in the SAURON sample. In contrast, 3:1 mergers do 
not produce any galaxy with a twist of more than 20° , be it 
with or without gas. In this respect 3:1 mergers are in better 
agreement with observations. 

Counter-Rotating Cores: We define CRCs as a subgroup 
of KDCs with a sudden kinematic twist of 180°. We find 
no CRCs in 3:1 remnants with or without gas and we find 
counter-rotating populations in two 1:1 coUisionless rem- 
nant, but which are not compact enough to pass as a core. 
For the stellar component of the 1:1 simulations with gas, 
we find at least four candidates for CRCs from the maps: 
11S4, 11S6, 11S1 3 and 11S14 and one in the gas 11GS16. 
iMcDermid et al.1 12006) found in the SAURON galaxies two 
dinstinct types of KDCs: slow rotators harbour large KDCs 
(kpc scale) with ages and metallicities similar to the main 
body of the galaxy and fast rotators have smaller (below kpc 
scale) and younger KDCs. Our simulations imply that stellar 
CRCs in slow-rotators are formed in 1:1 mergers when gas 
is present and, consistent with observations, as these CRCs 
consist of 'old' star particles in the simulation, they would 
not be detected as a distinct population, e.g. the CRC ap- 
pears old. Additionally the main body would rotate slower, 
if he has formed from a 1:1 merger. CRCs formed in the sim- 
ulations in the gas are much smaller and younger, again in 



agreement with observations, but we can hardly form them 
in 3:1 mergers. Most KDCs in fast rotators are actually 
CRCs, which might be a bias i n the observations (see discus- 
sion m iMcDermid et all 120061 ). still they seem too common 
place and merger simulations need to answer how to form 
them and maintain a high rotation in the main body. 
The V — hi Correlation: CoUisionless simulations show 
that there is generally a positive correlation between these 
quantities in the center of the remnants (NJB06). The pres- 
ence of gas would inhibit the population of box orbits, the 
orbit class responsible for positive correlations (see discus- 
sion in Sec. l4.3p . Contrary to what is seen in the coUisionless 
remnants, we see no correlation oiv — h^ in the center of any 
galaxy in the SAURON sample. In the outer parts there are 
a few examples with positive correlations: NGC 1023, NCC 
2699, NGC 4270 and particularly NGC 4526. This could b e 
an indication for a presence of bars (|Chung fc Bureaull2004l ). 
This argues (see NJB06) against pure coUisionless mergers as 
possible formation mechanism of rotating SAURON galax- 
ies, in agreement with the discussion on CRCs, see above. 
Opening Angle of Iso- Velocity Contours: In Sec. l4.5l we 
showed that on average 3:1 merger remnants with open iso- 
velocity contours have had less dissipation. This is because 
dissipation causes fiat discs and an overall more fiattened 
and axisymmetric potential. As more material now rotates 
in the equatorial plane, the iso-velocity contours are closing 
resulting in a smaller q. This result can not be easily com- 
pared in a statistical way to observations. Therefore we have 
chosen a pair of galaxies, which in the range of the SAURON 
observations have the same ellipticity e and the same v/a: 
these are NGC 7457 and NGC 4660, with e = 0.44 and 
v/a = 0.6 (|Cappellari et al.ll2005l ). Although these galaxies 
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appear in the v/a-e dia gram at the same posi tion, we find 
that the velocity maps of lEmsellem et ahl (|2004l ) of these two 
galaxies look very different: the velocity contours in NGC 
7457 are much more open than those for NGC 4660. This 
indicates that much more dissipation was present in the for- 
mation of NGC 4660, probably causing a flattened disc-like 
component in this galaxy. The ellipticity in the outer parts, 
from the RC3 (de Vaucouleurs, 1993), also is very different: 
0.41 for NGC 7457 and only 0.21 for NGC 4660. This gives 
support for the presence of a second flattened component in 
the middle of NGC 4660. For NGC 7457 we can not infer its 
intrinsic flattening, but it is likely to be c onstant with radius, 
since the ellipticity profile is constant (jPeletier fc Balcellsl 
1 19971 ). 

Polar Rings: As laid out in Sec. 14.71 we find one very pro- 
nounced polar ring in the merger 11GS9. This polar ring is 
present in simulations with or without gas. In the SAURON 
survey the only two examples of polar rings are NGC 2685 
and NGC 2768. B oth polar rings are only visible in the gas 
l|Sarzi et al.l l2006f ) . The radial extent of the polar ring in 
the merger is very similar to what is found in these two 
observations. There appears to be a fundamental difference 
between simulations and observations, as no stellar polar 
ring is found in the data. The statistics, however, is very 
small, so no strong conclusions can be drawn here. There are 
obser vational cases known of stellar polar rings (e.g. NGC 
4365 ISurma fc Benderl IT995I : iDavies et al.ll200ll ). but these 
are generally small, not going as far out as one effective ra- 
dius. 

Low a Rings: In some equal-mass merger remnants with 
gas we flnd ring- like depressions in a, mainly for face-on 
projections. Without gas we do not find them at all. Such 
a feature has been reported for one SAURON elliptical 
galaxy only: NGC 5813. I n the spiral galaxy NGC 4314 
l|Falc6n-Barroso et al.|[200^ ) a similar situation is visible: a 
star formation ring with a low gas velocity dispersion and a 
high stell ar dispersion. In a re cent survey of 6 nearby Seyfert 
galaxies (|Barbosa et al.ll2006l ) with IFU-GMOS, three of the 
galaxies, NGC2273, NGC3227 and NGC4593, show ring-like 
depressions in a in varying strength. If such low o rings are 
signatures of previous merging, it is interesting to note that 
most of them are found in active galaxies. 
Central a Drops: In almost all merger remnants which 
were formed from progenitors with a dissipational compo- 
nent, a significant fraction of the gas falls to the center, 
showing a very low aios in the maps. In contrast, the veloc- 
ity dispersion of the stars is increasing towards the center. 
In general the SAURON observations for early-type galaxies 
also show peaks in g in the center, in agreement with the 
simulations. There are a few objects with a central minimum 
in the velocity dispersion: these are NGC 4382 and NGC 
2768. For spiral ga laxies, recently ICanda et al. I l|2006l ) and 
iFalcon-Barroso et a l. (2006) found that dips in the central 
velocity dispersion are common, and that their frequency 
goes up with increasing Hubble type. They conclude that 
these dips are caused by central discs. In the case of the el- 
lipticals the stars nee d about IGy r to dynamically heat up 
after formation ( Wozniak et al.|[2003 ') . So we suspect that in 
these two galaxies the velocity dispersion dips are caused by 
young stars. For NGC 4382 this co uld be consistent with th e 
line-strength maps measured by (|Kuntschner et al] |2006| ) , 
but this idea does not work for NGC 2768 for which old 



populations are inferred in that paper. In this galaxy there 
must have been a different mechanism to keep the stars in 
the central region cold. This would be consistent with the 
center of some Sa galaxies, where velocity dispersion dips 
are found which are old (Peletier et al. in prep). 
Other a features: We found one example for a high a 
disk-like (llCl) structure caused by counter-rotation in 
the merger remnants. This feature is also rare in galaxies, 
NGC4473 being the only example in the SAURON survey 
(NGC4550 might also be comparable, but has a more over- 
all complex structure). A merger origin seems very plausi- 
ble for such a system although llCl has no net rotation and 
NGC4473 has still some. But the set-up of the merging discs 
gave them identical amount of rotation with opposite sign 
which is very unlikely to happen in nature, such that some 
net rotation might well survive a similar merging geome- 
try. NGC2549 and especially NGC3384 have dumbbeU struc- 
tures in their velocity dispersion maps and NGC3377 has a 
strong low g disk. All these galaxies are fast-rotating, which 
fits well with our results, that we find such disks predom- 
inantly in 3:1 mergers. However, we find them more abun- 
dantly than seen in EM0 4 The merger or igin of NGC2685 
has been proposed before (|Hau et al.lll999l ). We add circum- 
stantial evidence to this picture on the basis of the forma- 
tion of a double-peaked velocity dispersion structure in some 
equal-mass mergers. The dissipational component and the 
violence of the merging seem to be of the essence to puff-up 
stars and gas simultaneously. 

In general, coUisionless remnants show bigger opening 
angles of the iso-velocity contours, more boxy isophotes, ve- 
locity dispersion profiles without central peaks, no strong 
anti-correlation between v — h^. Remnants with a dissipa- 
tional component show smaller opening angles of the iso- 
velocity contours for the same ellipticities (inclinations), 
strong dips in the gaseous velocity dispersion profile and an 
anti-correlation between 11-/13. Additionally CRCs are more 
common in equal-mass mergers than in mergers of unequal 
mass. We conclude that globally the non-rotating subset of 
the representative SAURON sample of local galaxies agrees 
better with our 1:1 merger simulations, while the rotating 
subset can be reproduced by the dissipational 3:1 merger 
remnants. 



6 SUMMARY AND CONCLUSIONS 

We presented 2D maps of various moments of the LOSVDs 
of a large sample of 1:1 and 3:1 disc merger remnants with 
and without gas. Every remnant was resimulated with a dis- 
sipational component containing 10% of the luminous mass, 
allowing us to assess the influence of gas on the 2D fields. 
Additionally we performed a kinemetric analysis using the 
method devised by K06 to quantify properties such as kine- 
matic position angle or deviations from regular rotation. 

The difference between equal-mass and an unequal- 
mass merger remnants is no t only seen in pho tometric and 
global kinematic properties (Naab ot al. 199a), but also in 
the 2D kinematics (Bendo & Barnes 2000), which exhibits 
very different features. 1:1 mergers can lead to a merger rem- 
nant with low rotation. Also, orbit classes rotating around 
the major axis of the remnant can only be populated in 
significant numbers in 1:1 mergers, leading to kinematic 
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misalignment and twists, while 3:1 remnants rotate much 
faster and have almost no kinematic twists. All rotating 
coUisionless mergers fail to reproduce the observed v ~ h-j, 
anti-correlation, because the population of box orbits in the 
center is not inhibited (NJB06). 

The presence of a dissipative component of only 10% 
of the luminous disc mass in the progenitors can lead to 
a considerable change in the properties of 1:1 remnants. 
They are more round, kinematic misalignment is reduced, 
counter-rotating cores are forming and kinematic twist are 
more pronounced. Rare but observed features like polar gas 
rings form also in equal mass mergers. The effect of gas on 
3:1 remnants is less dramatic. KTs, CRCs and polar rings 
are not formed or at least must be very rare. A noticeable 
change is that the opening angles of the iso-velocity con- 
tours become smaller. As shown in NJB06, the v — hz anti- 
correlation is now much better reproduced, especially if the 
gas component is included in the analysis. 

Merger remnants show a variety of velocity disper- 
sion features. Significant central velocity dispersion dips are 
caused by infalling gas. The central star component in the 
gaseous runs, however, is heated. Edge on gas discs can ex- 
hibit lower velocity dispersion (as found in 3:1 remnants) or 
edge-on gas rings can produce double-peaked velocity dis- 
persion maps (as found in 1:1 remnants). Gaseous structures 
can imprint a ring-like depression in a even on the stellar 
component, which is also visible in hi. 

2D analysis of LOSVDs of N-Body simulations of galaxy 
formation are an additional tool to connect the formation 
history of galaxies with their observable kinematic fine struc- 
ture. We defer a full statistical analysis to future work, with 
a larger sample of galaxies mergers, which will also take star 
formation into account. 

The simulations and analysis presented her are a 
step forward in understanding the possible formation of 
ellipticals by mergers of discs. We focus on the influ- 
ence of a small dissipative component, which is impor- 
tant to understand its impact on the dynamical proper- 
ties of the remnants. However, additional physical pro- 
cesses like star-formation and feedback will change the 
detaile d properties especially at the center of the rem- 
nants jMihos fc Hernguistl ll996llSpringel et al.l 120051 ) and 
definitely have to be considered for simulations with higher 
gas f ractions (|Robertson et al.l l2004llSpringel fc Hernguistl 
I2OO5I ). We believe that the global effects of dissipation even 
with star formation will be very similar to the results pre- 
sented here. However, larger gas fractions might result in 
more massive stellar discs in the remnant. From a kinematic 
point of view disc merger remnants appear very similar to 
observed ellipticals, although questions regarding the age 
and metalicity of the stellar populations have to be ad- 
dressed in the future. A nalytical models of disc formation 
Naab fc Ostrikei|l2006l ) s well as semi-analytical modeling 



Khochfar fc Silk 2005) in combination with merger simula- 



tions l|Naab et al. .2006 ) can be used to place further con- 
straints on the disc merger hypothesis. 
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APPENDIX A: 2D VELOCITY MAPS AND 
KINEMETRIC DATA OF MERGER REMNANTS 

This is a selection of 2D velocity maps of aU the remnants. 
We have to restrict ourselves to the XZ projection, which is 
best suited for a kinemetric analysis, as it shows the highest 
rotational amplitude. For the equal-mass mergers we ana- 
lyze the maps by extracting the rotation curve along cir- 
cles, as the kinematic twists are too strong. Although we 
do not study projection effects, we can see the merger to 
merger variance of the kinemetric parameters, such as the 
kinematic misalignment in 1:1 mergers. A full list of the or- 
bital parameters for all merging geometries is given in Tab. 
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Figure Al. Top: Velocity maps for the coUisionless 1:1 remnants for the XZ projection. Bottom: Corresponding kinematic position 
angles extracted along circles. By definition kinematic PAs of 90° and 270° signify rotation which is aligned with the photometric major 
axis (indicated by straight black lines), correspondingly 0° and 180° indicate maximally misaligned rotation (red dot-dashed lines). 
Strong kinematic misalignment is present in almost all remnants. 
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Figure A2. Top: Velocity maps for the stellar component of the 1:1 remnants which had progenitors with a gas component. Bottom: 
Corresponding kinematic position angles. Definition for horizontal lines as before. Sudden kinematic twists give indications for kinematic 
decoupled components in the center of the remnants. In the outer parts the remnants are more aligned than their coUisionless counterparts 
(see Fig. lAll l 
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Table Al. FuUe List of merging geometries. For unequal-mass mergers the first number indicates the orientation of the more massive 
galaxy as i\ and u}\, the second number indicates the orientation of the more massive galaxy as i2 and ui2- 
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Figure A3. Top: Velocity maps of the same 1:1 remnants like in the previous figure, but with the gas particles converted into stars. 
Bottom: Kinemetric position angles as before. Overall misalignment is reduced, some KDCs are not visible anymore, as the gas dominates 
the center. 
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Figure A4. The fci term (bulk rotation) and the k^/k\ (deviation from regular rotation) are plotted with a commen x-axis. The 
low-rotation of 11C5, due to a high box orbit fraction in this remnant causes a high fcs/fci ratio. 

will be available under I http: / /www. usm.uni-muenchen.de] 

Figure A5. The ki (bulk rotation) profiles and the k^/ki ratio for the stellar component of remnants with gaseous progenitors. Some 
CRCs are clearly visible in 11S6, 11S13, 11S14 and 11S16 in both a sudden peak in k^/ki and in the double hump structure of ki. 
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Figure A6. Velocity maps for the coUisionless 3:1 remnants for the XZ projection. 

will be available under pittp: / /www. usm.uni-muenchen.de] 
Figure A7. Velocity maps for the stellar component of the 3:1 remnants which had progenitors with a gas component 
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Figure A8. Velocity maps of the same 3:1 remnants like in the previous figure but with the gas particles converted into stars 



